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We report a study of the decays B Dsi{2536)^ D^*\ where I)'*' is I)", D~ or D*~ , using a 
sample of 657 x 10® BB pairs collected at the T(45) resonance with the Belle detector at the KEKB 
asymmetric-energy e'^e~ collider. The branching fractions of the decays B'^ Dsi{2536)^ , 
I>si(2536) + L>- and B° D^i (2536) + I>*- multiplied by that of Dsi(2536)+ {D*°K+ + 
D*+K°) are found to be (3.97 ± 0.85 ± 0.56) x 10"'', (2.75 ± 0.62 ± 0.36) x 10"** and (5.01 ± 1.21 ± 
0.70) X 10"**, respectively. The ratio B{Dsi D'°K+)/B{Dsi D*+K°) is measured to be 
0.88 ±0.24 ±0.08. 



PACS numbers: 13.25.Hw, 14.40.Lb, 14.40.Nd 



The _Dsi(2536)+ is a narrow P-wave resonance for 
which the = 1+ assignment is strongly favored 1|. 
Although the Dsi{2536)^ was first observed in 1989 3|) 
its properties are still not well measured. In this anal- 
ysis we study the production of Dsi(2536)"*" in doubly 
charmed B meson decays, B — >■ Dsi{2b3Q)^ D''*\ where 
is either aD°, D' or D*' ^. The branching frac- 
tion measurements of the decays B — _Dsi(2536)+£'^*-' 
together with those of i? — >■ Z)(*^D^^^^ decays provide im- 
portant information to check the molecular hypothesis for 
the i:>*o(2317) and L'si(2460) partjcles First obser- 

vations of the B — > Z)si(2536)+D'^*^ decay modes have 
been reported by BaBar 0, 01 ■ An upper limit on the de- 
cay -J> Dsi{2536)+D*- was also obtained by Belle 
which is consistent with the BaBar measurement. 

This analysis is based on 605 fb~^ of data collected at 
the T{4S) resonance with the Belle detector at the 
KEKB asymmetric-energy e~^e~ collider [loj . which cor- 
responds to 657 X 10^ BB pairs. The Belle detector is a 
general-purpose spectrometer with a 1.5 T magnetic field 
provided by a superconducting solenoid. A silicon vertex 
detector and a 50-layer central drift chamber are used to 
measure the momenta of charged particles. Photons are 
detected in an electromagnetic calorimeter consisting of 
CsI(Tl) crystals. Particle identification likelihoods Ck 
and are derived from information provided by an ar- 
ray of time-of-flight counters, an array of silica aerogel 



Cherenkov threshold counters and dE/dx measurements 
in the central drift chamber. Two inner detector config- 
urations were used. A 2.0 cm radius beampipe and a 3- 
layer silicon vertex detector were used for the first sample 
of 152 X 10^ BB pairs, while a 1.5 cm radius beampipe, a 
4-layer silicon detector and a small-cell inner drift cham- 
ber were used to record the remaining 505 x 10^ BB 
pairs. 

All charged tracks are required to have a distance of 
closest approach to the interaction point (IP) in the plane 
perpendicular to the beam axis smaller than 2 cm and 
smaller than 5 cm along the beam axis. Charged kaon 
and pion candidates are required to be positively identi- 
fied. The Kg candidates are reconstructed in the tt+tt" 
mode with the requirement |M^^ — m^^oj < 15 MeV/c^ 
(3cr), where rrij^o is the Kg mass pj]. Requirements on 

the Kg vertex displacement from the IP and on the dif- 
ference between the vertex and Kg flight directions are 
applied. No pion identification is required for the pions 
from Kg candidates. A mass- and vertex-constrained 
fit is applied to improve the four-momentum measure- 
ments of Kg candidates. Photons are reconstructed in 
the electromagnetic calorimeter from showers that are 
not associated with charged tracks with energies larger 
than 50 MeV. Combinations of two photons are consid- 
ered to be tt'^ candidates if their invariant mass lies within 
±15 MeV/c^ (3(t) of the 7r° mass [ll|. To improve their 
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momentum resolution, all tt" candidates are fitted with a 
TT^ mass constraint. Continuum e+e^ — >■ qq backgrounds 
{q = u, d, s, c) are suppressed by requiring the ratio of 
the second and zeroth Fox- Wolfram moments [ll| to be 
smaller than 0.3. 

Candidate D°'s are reconstructed using five decay 
modes: K~7r~^, K'tt'^tt^tt' , K~Tr~^7T'^, Kgir^ir' and 
K'^K". The is reconstructed via its decay into 
if+TT^TT^. The selected combinations are constrained to 
a common vertex and the /n.d.f. of the vertex fit is re- 
quired to be smaller than 25. A ±15 MeV/c^ {2a) mass 
window around the mass [l[ is used to select can- 
didates for the £>° K-tt+t:^ mode and ±12 MeV/c^ 
(« 3cr) for the other D decay modes. To reduce the large 
combinatorial background in the — > if ~7r+7r° decay 
mode, the energy of each photon from tt*^ is required to 
be greater than 100 MeV. For the D^^s coming directly 
from B decay, only the cleanest modes, if+7r~, if^Tr+TT" 
and K^K~, are used. A mass- and vertex-constrained 
fit is applied for all D candidates. The _D*+ is recon- 
structed via its decay to D^t^sIow "^'^ improve the tt^q^ 
momentum resolution, and thus the mass resolu- 

tion, the TT^o^ is constrained to the decay vertex, 
which is obtained by fitting the to the IP profile. The 
£)*+ candidates with invariant mass within ±2 MeV/c^ 
(4ct) of the D*+ Q mass are selected. The D is re- 
constructed using the D°7r° and decay modes with 
invariant mass windows of ±3 MeV/c^ and ±12 MeV/c^ 
{2a), respectively. In both cases a mass-constrained fit 
is performed for D*'^ candidates. 

The £)si(2536)+ meson is reconstructed in its domi- 
nant decay modes: D*^K'^ and D*^K^. The invariant 
mass of the -Dsi(2536)^' candidates is required to be less 
than 2.58 GeV/c^. Combinations of Dsi(2536)+ and a 
second charm meson {D^*^ = D°, D~ , D*~) with 
opposite charm flavor are considered as B meson candi- 
dates. B candidates are identified using the energy dif- 
ference AE — Eb — -Ebcam and beam-energy constrained 
mass Mbc = V-^^boam ' where Eb and pg are the B 
candidate energy and momentum in the center-of-mass 
system, respectively. We require |Ai<^| < 0.1 GeV and 
Mbc > 5.22 GeV/c^ to preselect B candidates. After 
these selections the average candidate multiplicity per 
event is 1.6. A single candidate per event is chosen based 
on the smallest of the Z)'^*^ candidates mass devia- 
tions from their nominal values. The AE versus Afbc 
scatter plot for the candidates with Dsi{2536)'^ mass 
within ±5 MeV/c^ {3a) of the nominal value [l| is pre- 
sented in Fig. [T] left, for the sum of all studied decay 
modes. A signal box is defined as \AE\ < 0.02 GeV and 
Mbc > 5.27 GeV/c2 (region "a"). The AE - Mbc two- 
dimensional sideband is defined as Afbc < 5.26 GeV/c^ 
or \AE\ > 0.03 GeV (region "b"). The Mbc {AE) pro- 
jection is shown for the events with |A£'| < 0.02 GeV 
(Mbc > 5.27 GeV/c^) in Fig. Reenter (right). 

A two-dimensional binned likelihood fit is performed 
to the AE — Afbc distribution and yields 151 ± 15 signal 
events. The signal yield obtained from the AE — Afbc 



fit contains a contribution from non-resonant B — >■ 
D*KD^*^ decays that have the same final state parti- 
cles as the signal. Therefore to extract the yields of 
each decay mode, we perform an additional binned like- 
lihood fit to the £'si(2536)+ mass distributions for the 
events in the AE — Afbc signal box (Fig. [2]). In total, 
there are nine reconstructed decay modes corresponding 
to three B meson decay modes, i3+ f)si(2536)+-D°, 
^ f:>^i(2536)+f?- and B° ^ f?,i(2536)+L>*-, and 
three f?si(2536) meson decay modes, f?si(2536)+ 
D*°{D°j)K+, D*°{D°TT°)K+ and D*+ {D°tt+)K°. The 
latter are related by the known branching ratios R' = 
B{D*° D°n°)/B{D*" ^ D°j) = 1.74±0.13 and f?" = 
B{Dsi{2536)+ -> D*°K+)/B{Dsi{2536)+ ^ D*+K°) = 
1.36 ± 0.20 [l|. The shapes of the signal distributions 
are taken to be a non-relativistic Breit-Wigner functions 
convolved with the f?si(2536)+ mass resolution func- 
tion. The Breit-Wigner mass and width of f?si(2536)"'" 
are floating parameters common to all modes. The 
f3si(2536)"'" mass resolution is parametrized by a double 
Gaussian; the resolution parameters and reconstruction 
efficiencies are obtained from Monte Carlo (MC) sim- 
ulation and are summarized in Table HI The shape of 

TABLE I: The products of the total efficiencies and inter- 
mediate branching fractions, eB, and the Dai (2536)^ mass 
resolution parameters, where amain (ftaii) is the width of the 
narrow (wide) Gaussian component, and /tail is the fraction 
of the wide component. 



Dsi mode: 


D*"{D"'y)K+ 


f3*"(D'V'^)A-+ 


D'+Ks 


eB{DsiD"), 10-* 


1.45 


1.18 


0.74 


eB{D,^D-), 10-* 


1.97 


1.62 


1.00 


eB{DsiD'-), 10-* 


0.93 


0.75 


0.32 


(Tmain, McV/c^ 

(Ttaii, MeV/c^ 


1.02 


1.01 


0.95 


5.44 


3.61 


2.43 


/tail 


0.23 


0.34 


0.15 



the background is parametrized with an ARGUS func- 
tion jT5] and fixed from a fit to events in the AE — Afbc 
sideband. This background function includes both the 
non-resonant component and the combinatorial back- 
ground since their shapes are expected to be similar. All 
nine distributions are fitted simultaneously. The ratios 
of signal yields in different f?si(2536)"'" decay modes are 
fixed using their relative branching fractions, R' and f?", 
and MC reconstruction efficiencies. The number of signal 
events for each of the B decay modes is fioated in the fit as 
a sum of the signal events in all Z)si(2536)+ decay modes; 
the background normalization is free for each distribu- 
tion. The results of the fit are presented in Table HIl The 
significance is calculated from S — \/ —2 ln(£o/'Cmax), 
where and £max are the maximized likelihoods with 
the signal yield fixed at zero and left free, respectively. 
The L»^i(2536)+ mass is found to be 2534.1±0.6 MeV/c^ 
and the width is F = 0.75±0.23 MeV/c^ consistent with 
the current world-average values in p|. To provide a 
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FIG. 1: The scatter plot of /S.E vs Mbc (left), where (a) is the signal region and (b) is the sideband, and projections on 
A_E for events with Mbc > 5.27 GeV/c^ (center) and Mbc for events with |Ai?| < 0.02 GeV (right). The Mbc distribution is 
parametrized with a single Gaussian for the signal and an ARGUS function flU for the background. The Ai? is parametrized 
with a double Gaussian and a linear polynomial for the signal and background, respectively. The results of the fits are shown 
by the superimposed curves. 



measurement of the £'si(2536)+ branching-fractions ra- 
tio another fit is performed with a floating R" parameter. 
From the fit we obtain R" = 0.88 ± 0.24, consistent with 
other measurements 

TABLE 11: Summary of the fit results: event yields, branching 
fractions, and statistical significances. 



B decay mode 


N 


B, IQ-^ 


S 


Dsi(2536) + (D*Jf)_D" 
Dsi(2536) + (D*A')_D- 
Dsi(2536) + (I>*A')D*- 


42.5 ± 8.7 

40.2 ± 8.6 

33.3 ± 7.6 


3.97 ± 0.85 
2.75 ± 0.62 
5.01 ± 1.21 


7.0cr 
6.9cr 
6.3cr 



The main systematics for the branching fraction mea- 
surements is from the tracking efficiency. A 1% system- 
atic error is assigned for each track and an additional 3% 
error for each low momentum track, which are summed 
linearly. For the kaon identification a 1% systematic un- 
certainty is assigned for each kaon track. The contribu- 
tions of the systematic uncertainties of the and 7/7r° 
reconstruction efficiencies to the overall systematics are 
estimated to be 1% and 3%, respectively. The width of 
the narrow Gaussian component of the signal mass reso- 
lution is increased by 20% to obtain the systematics due 
to the poorer resolution in data compared to MC simu- 
lation. To obtain the systematics due to the imperfect 
background shape description, we vary the shape of the 
function describing the background: the parameters of 
the ARGUS function are varied within their errors; we 
also fix the shape to the one obtained from the generic 
MC simulation, or use a square root function instead of 
the ARGUS function. The largest difference in the re- 
sults is treated as the systematic uncertainty in the back- 
ground shape description. A possible contribution from 
inclusive B Dsi{2bm)+X and qq Dsi{2bm)+X 



production is checked by examining the Z?si(2536)^ mass 
distribution in a BB MC simulation from which the sig- 
nal has been removed, and in data in the AE' — Mbc 
sideband; no signal is found in both cases, so an upper 
limit on the yield is used as the systematic uncertainty 
due to the peaking background. Another systematic un- 
certainty is from the errors in the fractions R' and R" of 
the Z?*° and Z?si(2536)+ subdecay modes, respectively. 
These ratios are varied in the fit within their errors, and 
the difference in the fit results are assigned as a system- 
atic uncertainty. All the individual systematic errors, 
shown in Table IIIIl are summed in quadrature. Since 
the total systematic uncertainty is nearly symmetric, the 
maximum of the positive and negative errors is taken as 
the final systematic uncertainty. The errors in the results 
are dominated by statistical uncertainties. 



TABLE III: Relative systematic errors for the branching frac- 
tions and the ratio of the branching fractions, B{Dsi 
D''°K+)/B{Dsi ^ D*+K°), in %. 



Source 


B{DsiD") B{DsiD- 


) BiDsiD- 


R" 


Tracking 


±9 


±9 


±9 


±3 


Particle ID 


±3 


±3 


±3 


±1 


7/7r0 


±3 


±3 


±3 


±3 


Ks 


±1 


±1 


±1 


±1 


MC resolution 


±3 


±3 


±3 


±6 


BG shape 


±4 


±4 


±4 


±3 


Peaking background 


+0 
-3 


+0 
-3 


+0 
-3 


±3 


B(D(*') 


±7 


±5 


±7 


±2 


R' 


±1 


±1 


±1 


±1 


R" 


±1 


±1 


±1 




N{BB) 


±1.5 


±1.5 


±1.5 




Total 


±14 


±13 


±14 


±9 
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FIG. 2: Dsi(2536)+ mass distributions for: a), b), c) B+ L>si(2536)+73°; d), e), f) B° D si (2536) + D' and g), h), i) 
B° Dsi{2536)+ D*- final states, foUowed by L>si(2536)+ decays to a), d), g) D*°{D°'y)K+; b), e), h) D*°{D'-\°)K+ and c), 
f), i) D*+{D'\+)K°g. The points with error bars are the data, while the curves show the fit result. 



In summary, we report a measurement of the branch- 
ing frac_tions for the decays B -J> Dsi(2536)+Z)(*\ 
where D'-*'^ is D°, D- or D*- . From a simultane- 
ous fit to all B and I?si(2536)+ decay channels we 
measure B{B+ L>^i(2536)+L'°) x i3(i:>^i(2536)+ ^ 
{D*^K+ + D*+K^)) = (3.97 ± 0.85 ± 0.56) x 
10-^ B{B^ Dsi{2536)+D-) x S(i:»si(2536)+ ^ 
{D*°K+ + D*+K^)) = (2.75 ± 0.62 ± 0.36) x lO"'^ and 
B{B° Dsi{2^iQ)+D*-)xB{Dsi{2bm)+ {D*^K+ + 
D*+K°)) = (5.01 ± 1.21 ± 0.70) X lO"''. The ratio 
B{Dsi D*°K+)/B{Dsi D*+K°) is measured to 
be 0.88 ±0.24 ±0.08. The first error is statistical and the 
second one is systematic. The obtained results are con- 
sistent within errors with the previous measurements 0] . 

Using the latest measurements of the B — > D'^*^D^*X. 
branching fractions pj we calculate the ratios discussed 
in Q: 



Rdo 



Rd*o — 



Rdi — 



Rd' 



BjB ^ DD*,i2317)) 

B{B DD,) 
B{B D*D*„(2S17)) 

BiB D*D,) 
B{B DDsi{2460)) 

B{B ->■ DDI) 
B{B (2460)) 
B{B D*D*) 



0.10 ±0.03, 



= 0.15 ±0.06, 



= 0.44 ±0.11, 



0.58 ±0.12. 



In addition, the same ratios are calculated for B 
D<-*^D si{2536)+ decays using combined BaBar [71 and 
current results: 



Rdi 



R 



D-V 



BiB DDsii2536)) 

B{B DDI) 
B{B D*Dsii2536)) 
B{B D*D*) 



0.049 ±0.010, 



= 0.044 ±0.010. 



In these calculations it is assumed that the decay modes 
£1*0(2317)+ ^ D+TfO and 15,1(2536)+ {D*°K+ + 
]j*+X^) are dominant. 

According to [1, [1] , within the factorization model and 
in the heavy quark limit, these ratios should be of order 
unity for the L»*o(2317) and Dsi(2460), whereas for the 
Z?si(2536) they can be very small. From the above ra- 
tios we can conclude that while the decay pattern of the 
£'si(2536) follows the expectations, the new Dgj states 
are either not canonical cs mesons, or this approach does 
not work for these particles. 
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